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de Haas-van Alphen oscillations in quasi-two-dimensional underdoped cuprate 
superconductors in the canonical ensemble 
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We calculate the de Haas- van Alphen (dHvA) effect waveform using the canonical ensemble for dif- 
ferent Fermi surface scenarios applicable to the underdoped cuprate superconductor YBa2Cu306.5, 
in which quantum oscillations have recently been observed. The harmonic content of the dHvA 
waveform of the principal Fa ~ 500 T frequency is consistent with the existence of a second thermo- 
dynamically dominant section of Fermi surface that acts primarily as a charge reservoir. Oscillations 
in the charge density to and from this reservoir are shown to potentially contribute to the observed 
large quantum oscillations in the Hall resistance. 

PACS numbers: 71.45.Lr, 71.20.Ps, 71.18.+y 
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Layered electronic structures are known to be favorable 
for unconventional superconductivity in several families 
of compounds [l|, l3j |3| ■ The quasi- two-dimensional Fermi 
surfaces of layered systems also yield magnetic quantum 
oscillations notably different from those in conventional 
three-dimensional metals. When the interlayer hopping 
energy of the quasiparticles t^ becomes less than the cy- 
clotron energy heB/ni*, increased Landau level degener- 
acy causes the quantum oscillations to depart from the 
Lifshitz-Kosevich (LK) theoretical form 4]. The extent 
to which the chemical potential oscillates plays a cen- 
tral role in determining the de Haas- van Alphen (dHvA) 
waveform [a, 0, Q- In effect, chemical potential oscilla- 
tions become a valuable tool for probing the thermody- 
namic mass of additional Fermi surface sections [6J, 01 
through their effect as a 'charge reservoir,' and have 
been shown to play an important role in layered or- 
ganic metals [7|, |9|, [lO|. In this paper, we propose two- 
dimensionality to play an equally important role in deter- 
mining the dHvA waveform and Hall resistance oscilla- 
tions in underdoped cuprate superconductors |lll.ll2l.ll3l|. 
Importantly, the extent to which the chemical potential 
oscillates can help distinguish between Fermi surface sce- 
narios consisting of an isolated pocket or of multiple sec- 
tions. 

The motivation for this study stems from the recent 
observation of a possible harmonic and second dHvA fre- 
quency in YBa2Cu306.5i [111. These features are dis- 
cernible in Fig. [1] on performing a Fourier transform 
after subtracting a damped sinusoidal fit to the domi- 
nant Fa '^ 500 T component [ll|. The lower frequency 
in Fig. [T] is twice Fa, which is consistent with it be- 
ing the harmonic i^2a- Furthermore, on fitting the en- 
tire dHvA waveform to a sum A — ^^ Ai sm{2TTFi/B + 
(pi) exp(— Fi/B) over frequencies i in Fig. [2^, the Fa and 
F2a oscillations can be seen cross through zero at sim- 
ilar values of the magnetic field— thus exhibiting the 
key distinguishing feature of the dHvA effect in a two- 
dimensional metal [14| . 

To model the dHvA oscillations in YBa2Cu306.5i pTj . 
we therefore consider simplified quasi-two-dimensional 
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FIG. 1: dHvA oscillations suffer diffraction when Fourier 
transformed over a finite interval in magnetic field, with the 
side lobes sometimes interfering destructively with less promi- 
nent features. This situation can be remedied by subtract- 
ing fits to prominent features, such as the Fa ~ 500 T fre- 
quency [12| ]. prior to Fourier analysis. The black curve shows 
the result of such a procedure, yielding frequencies F2a and 
Fp . The red curve shows the same region of the Fourier trans- 
form prior to subtraction while the cyan curve shows the resid- 
ual after subtracting fits to all three frequencies. 



dispersions e = r]ih (k^ + ky)/2m* + 2ic,i(l ~ cos[kzc]) 
for each Fermi surface section, where rji = -fl for 
electrons and —1 for holes, c is the interlayer spacing 
and tc,i <^ eB/m*. A large magnetic field along z 
causes the density-of-states for each section gi[e,B] = 
(rjieBTi/'K'^hc) XllTC^ilH^ + ^1)"^ to become a series of 
Landau tubes v where e^ = s—eo^i — r]iheB{i' — l/2)/m*± 
Asi [a,[3|- Here, Asi represents a possible splitting of the 
Landau levels caused by Zeeman splitting |15j . bilayer 
splitting or ic,i [16] . Scattering from random defects and 
impurities is assumed to broaden the Landau tubes into 
Lorcntzians of width 2Ti ~ h/ri |l3| . 

Since the model treats values of ic,i ^ eB/m*, the 
entire Fermi surface contributes to the dHvA effect, re- 
quiring the magnetization M = —dH-p/dB\is!.T to be cal- 
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FIG. 2: Decomposition of the experimental data (a) and nu- 
merical model simulations (b,c and d) into their constituent 
Fourier components. Note that the sign of the magnetization 
is opposite to that of the torque presented in Ref. [ll|l. The 
simulations are performed for a two-dimensional Fermi surface 
with background charge reservoir of different relative masses. 
The amplitudes of the harmonics and Fjs are increased x 5 
for clarity. 



culated numerically using the multiband canonical en- 
semble 0, [a, 0, Q- The free energy is given by Hp = 
T,zI-oo^9i[s,B]{l ~ exp[?7,(^ - e)/kBT)])-'^de, requir- 
ing the chemical potential fi to be determined by fixing 
the total number of particles N = J2i I- 5i[^i^](l ~ 
exp[r],{fi-e)/kBT)])-^de. 

We begin by considering a closed pocket (a) in the 
presence of a background uniform density-of-states ^rcs — 
m*^g/Trh c, which functions only as a charge reservoir 
(corresponding either to an open Fermi surface sec- 
tion or to a closed section with r^^ ^ ^a^)- Fig- 
ures [2}D-d show calculations of the magnetization for 
background reservoirs of different effective masses, us- 
ing Fa = 499 T to match the phase of the oscillations 
measured in Fig. [2^ (l8|.and m* = 1.9 me [il|, Ae^ = 
and T~^ — 2.5 x 10^^ s~^ |19f . The scenario depicted 
in Fig. [^b corresponds to that of an isolated pocket (or 
two or more symmetry related pockets) with no c harg e 
reservoir present (i.e. m*^^ = 0,see Refs. p^. [2(1 [2ll. [2a| ) . 
However, on comparing Figs. [2^ and b, the amplitude of 
the i<2a frequency is this scenario can be seen to have 
the opposite sign of that observed experimentally [23| . 
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FIG. 3: Colored contour plot of the harmonic ratio M2a/Ma 
in the simulations versus 771*^5 /m* and r~^ for the same aver- 
age field range (41 T) as in Fig. [2] The M2c«/Mq = -0.03 con- 
tour corresponds most closely to that representing the data in 
Fig. EK while the dotted lines represent the likely upper and 
lower bounds of the experimental scattering rate (see text). 



In contrast, an additional charge reservoir of thermody- 
namic mass m*j,g > m* (e.g. Fig. [2Ji) can lead to simu- 
lated F2a oscillations with the same sign as the experi- 
mental data. 

Figure [3] shows the harmonic ratio Mia I Ma for the 
same average field interval as in Fig. [2] plotted for a 
variety of t~^ and 771*^3/771^ values. The uncertainty 
in the contribution to the effective scattering rate from 
the superconducting pairing potential within the vortex 
state [2J| requires us to consider a range of t~^ val- 
ues. A precise estimate would require measurements 
at magnetic fields [IqH ^ p^^Uci ^ 60 T [ll|. Since 
the contribution to t~^ from the vortex state is mag- 
netic field-dependent [2J], estimates made from the gra- 
dient of a simple Dingle plot below H^i typically over 
estimate the total extent of the scattering [ll|. The 
valu e T~ ^ « 5 X 10^^ s^^ obtained from such an anal- 
ysis [111 [l3, [i3 therefore can only be regarded as an 
upper bound (rightmost dotted line in Fig. [H]). A lower 
bound estimate can be obtained from the extent to which 
the oscillations are damped between the end limits of 
the interval in magnetic field over which experiments are 
performed. From the oscillations of Jaudet et al. [1J| 
measured between ~ 30 and 60 T, we obtain a lower 
bound of T~^ « 2 X 10^^ s~^ (leftmost dotted hne in 
Fig. [21). The region bounded by the dotted lines and 
by the uncertainty in the experimental harmonic ratio of 
Mia/Ma = -0.03 ± 0.01 therefore provides us with the 
permissible range of parameters for modeling the dHvA 
effect, roughly corresponding to Ta^ ^ 3.5 x 10^^ s~^ 
and m*^^/m*a ^ 1.5. 

Based on the above comparisons, rather than exist- 
ing in isolation, the a pocket appears to coexist with 



a charge reservoir of thermodynamically greater mass. 
Open (or quasi-one-dimensional) Fermi surface sections 
originating from the Cu-0 chains or from Fermi surface 
reconstruction by a coUinear density-wave ordering with 
Q — (tt, (1 ± 2S)-k) J5j provide potential candidates for 
a uniform reservoir of the type necessary to generate a 
waveform Hke that in Fig. [5^. This would then require 
the Ffi oscillations to be the product of magnetic break- 
down tunneling [llj. 

An alternative possibility is that the reservoir is com- 
posed of closed Fermi surface pockets of significantly 
greater thermodynamic mass or scattering rate than the 
a pocket so as to appear as relatively weak features 
in the dHvA signal. With its heavier effective mass of 
m^ « 3.8 mo, the Ffi frequency in Fig. [2^ (reproduced 
in Fig. [4^) could correspond to this reservoir. A simu- 
lation including both a electron and /3 hole pockets in 
Fig. [4|3 with similar frequencies and effective masses to 
those reported in Ref. [llj] can produce absolute ampli- 
tudes in similar proportions to those in Fig. H^ on set- 
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= 2.1 X 10^ 



and 



3.1 X 10^ 



However, agreement with the sign of the observed /? fre- 
quency oscillations and the magnetic field dependence of 
its amplitude requires an additional effect, potentially 
provided by a fixed splitting term Ae^ = 0.84 meV 
we introduce into the simulation in Fig. [It (also using 
Tg^ = 2.6 X 10^^ ^~^)- This splitting could correspond 
to effects such as an interlayer tunneling tc.p = 0.42 or 
residual bilayer splitting [ll| . Such a value for Ae^j would 
give rise to a beat with nodes at ~ 18 and 55 T— presently 
outside the field range of the torque experiments. The 
reduction in the amplitude of the Fq oscillations with in- 
creasing field in YBa2Ba306.5o JLi] could be consistent 
with a node at higher magnetic fields. 

Another proposed scenario, corresponding to recon- 
struction of the Fermi surface by Q = (tt, tt), consists 
of hole pockets situated at k = (7r/2, 7r/2) of nearly dou- 
ble the fc-space area and frequency of proposed electron 
a pockets situated at k = (tt, 0) [2^. Since there are 
two hole pockets for every electron pocket in this sce- 
nario (therefore contributing twice as much to the dHvA 
signal and thermodynamic mass), the effective mass or 
scattering rate of the hole pocket would need to be sev- 
eral times that of the a pocket in order for its amplitude 
to be so much weaker in the experiment (e.g. requiring 
m* > 6 me or r"! > 7 x IQi^ g-i). 

Quantum oscillations in the Hall resistance [l2] could 
provide further clues as to the nature of the charge reser- 
voir. Doiron-Leyraud et al. [12] find quantum oscillations 
that are a proportionately larger fraction of R^y than 
Rxx, suggesting a possible dominant oscillatory contri- 
bution to the Hall resistivity from oscillations in the car- 
rier density N^ (and that of the reservoir) rather than 
axx or ayy [28|. Oscillations in the carrier density result 
from the oscillatory flow of charge back and forth between 
the pocket and the reservoir given by N^ — N^os, where 






J^^ga[e,B]{l - exp[7y„(^ - e)/kBT)]) Me and 
A'rcs.o — Na + Na and Ni refers to the non os- 
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FIG. 4: Data re-plotted from from Fig. [2^ (a) together with 
model calculations (b and c) for a Fermi surface consisting 
of a small electron pocket (a) and large hole pocket (/3). b 
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FIG. 5: (a) Oscillations in the Hall effect and magnetiza- 
tion from Refs. [13. [l3]. The data have been rescaled by an 



exponential factor e 



-7/S 



to facilitate comparison. (b)The 



Hall resistivity due to A^^ and A^rcs calculated using a two 
band model (assuming the a pocket to be electron-like) with 
non-oscillatory Drude conductivities determined by trans- 
port scattering rates as indicated. The black, red and blue 
curves have been renormalized by 42.2, 12.0 and 2.3 re- 
spectively. The inset shows a Fourier transform of Na for 
47 < ^oH < 62 T (roughly corresponding to the experimen- 
tal range [l2|) for a reservoir consisting of a closed pocket (/3) 
using the parameters as used in Fig.|4};. Consistent with Hall 
resistance experiments, the /3 frequency is not a prominent 
feature in Na- 



dilatory zero field values. Consistent with a dominant os- 
cillatory carrier density contribution |23| , the phase of the 
oscillations in the experimental Hall resistance [1^ and 
magnetic torque [l^ are aligned to within « 0.7 ^ 7r/2 
in Fig. [5^. 

Figure [5)3 shows the result of inserting the respec- 
tive Hall coefficients Ra = 1/eNa and i?rcs = 1/eA^rcs 
into a simple two-band expression for the Hall coefficient 
-Rh = (-RqCTq + -R/30'res)/(o'a + cTiog)^ On extending the 
simulations to higher magnetic fields. We neglect os- 
cillations of a^x = o'yy by assuming simple Drude ex- 
pressions axxA = eNi^oTi^t/fn*, where t~^ < t~ is the 
transport scattering rate. The oscillations are strongest 
(black curve) compared to the background on assuming 
a uniform transport scattering rate in the model includ- 
ing both a and f3 pockets from Fig. |4):, due to the very 
similar magnitudes of Racr^ and R,:esO'^cs ™ this case. 
They are also found to be prominent (red curve) on in- 
serting T~^ = T~^ from Fig.HJ:;. The oscillations become 
rather small (blue curve), however, if the reservoir has 
a poor conductivity by setting t"^ ^ > 7 x 10^^ s~^ or 
JTires ^ 6 rric as would be implied by the Q = (tt, tt) sce- 
nario discussed above, or if i?res = as for an ideal open 
Fermi surface. The oscillations vanish in the case of an 
isolated pocket with no reservoir (green curve) . It should 
be noted, however, that the present simulations cannot 



reproduce the reported magnitude of the Hall resistiv- 
ity [12] ■ This would require a single isolated a electron 
pocket with no Kramers degeneracy, the additional pres- 
ence of a very small much higher mobility electron pocket 
or a possible vortex liquid contribution [29|. 

In summary, by comparing the measured dHvA wave- 
form of YBa2Cu306.5i with canonical ensemble simula- 
tions, the sign of the harmonic of the dominant oscilla- 
tory Fa contribution suggests the coexistence of the a 
pocket with a thermodynamically dominant (i.e. heavier 
mass) charge reservoir. This finding is consistent with 
the the existence of multiple carrier sections [111 . |23 | . A 
simulation including both electron and hole closed pock- 
ets in the canonical ensemble appears to consistently ex- 
plain the main aspects of both the dHvA waveform and 
the phase of the oscillatory Hall resistance on allowing 
for small adjustments in the transport scattering rates. 
In-situ measurements of magnetization and Hall resistiv- 
ity should ultimately enable the relative contributions of 
the carrier density and diagonal conductivity to be better 
resolved. 
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